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We compute the Casimir interaction energy between two perfectly conducting, 
concentric cylinders, using the mode-by-mode summation technique. Then we 
compare it with the approximate results obtained using the proximity theorem 
and a semiclassical approximation based on classical periodic orbits. We show 



o 

o . 

^^H 1 that the proximity theorem with a particular choice for the effective area coincides 

(-H ' with the semiclassical approximation and reproduces the exact result far beyond 

r^ ' its expected range of validity. We also compute the force between slightly eccentric 

^_ cylinders and discuss the advantages of using a cylindrical geometry to measure 

the Casimir force. 



1 Introduction 



Many of the recent experiments that measured the Casimir force between con- 
ductors involved the interaction between a plane and a sphere,!^ with only one 
exception where the force between parallel plates was measuredP' The plane- 
>0 , sphere configuration is more appealing from an experimental point of view, 

mainly because the alignment problem that complicates the measurement for 
the parallel plates is not present. However, from a theoretical point of view, 
■~V ' the inconvenience is that there is no exact evaluation of the force between a 

"^ - sphere and a plane. The comparison between theory and experiment is done 

with the "proximity" approximation for the force. 
S . It is then of interest to explore other geometries for which the Casimir 

C^ ' force can be computed exactly. This would be useful both for testing the 

^ , validity of the approximations and, eventually, for proposing new experiments. 

S? ' In what follows we will be concerned with the Casimir interaction between 

two long, perfectly conducting cylindrical shells of radii a and b and length 
K/i L:^ a,b. 

5t 1 2 Exact Casimir energy 

In order to define properly the Casimir energy we introduce a cutoff a as 
follows 

Ec(.a) = ^ Y^ (e-^-'-w.-e-^^'-w,). (1) 

p—n.m,kz 

The exact Casimir energy Ec is the limit of Ecic) as cr —)■ 0. In our defini- 
tion we take Wp as the eigenfrequencies of the electromagnetic field satisfying 
perfect conductor boundary conditions at r — a,r — b and r = R, and Wp are 



those corresponding to the boundary conditions at r = i?i , r = i?2 and r — R, 
in the Hmit R > R2 > Ri ^ a > b. The parameters Ri , R2 and R define the 
reference vacuum. 

In cylindrical coordinates, the eigenfunctions are of the form 

h„k^ [t, r, 6, z) = e(-™-^*+^'=^-+^"'')i?„(Ar) , (2) 

where the function i?„ is a combination of Bessel functions satisfying the 
perfect conductor boundary conditions. These boundary conditions define the 
possible values of the constant A. The eigenfrequencies are w^fc^ = c\/k1 + A-^. 
The frequencies of the TE modes are defined by 

F^^{z,a,b,R) = Jn{za)[Jn{za)N„{zb) - J„(z&)A^„(za)] 

X [J^{zb)NnizR) - JnizR)Nnizb)] = . (3) 

The frequencies of the TM modes involve derivatives of the Bessel functions, 
and one can define an analogous function F^^^ . In the set of quantum num- 
bers p = {71,171, kz) appearing in Eq. Q, m denotes the different solutions 
Xnm of both i^J^(z, a, b,R)^0 and F^^'{z, a,b,R)^0 . 
From Cauchy's theorem it follows that 

J- f dzze-"^^\nf{z)=Y,x,e-^^^ , (4) 

2m Ju dz ^—' 

where /(z) is an analytic function within the closed contour C, with simple 
zeros at x\,X2, ■■■, within C. Wc use this result to replace the sum over m in 
Eq.(^) by a contour integral in the complex plane. "^ We findl^ 

(5) 

where 

pBM. ,, y F^^{z,a,b,R)F™{z,a,b,R) 

^" ^^'"'°^ Rur3.^ooFTE[z,Ri,R2,R)F™{z,R^,R2,R)- ^' 

An adequate contour for the integration in the complex plane is a circular 
segment Cr of radius T and two straight line segments forming an angle (f) 
and -K — (j) with respect to the imaginary axis. The nonzero angle is needed 
to show that the contribution of Cr vanishes in the limit F — > oo when tr > 0. 
It proves to be convenient to compute the difference between the energy of 
the system of two concentric cylinders and the energy of two isolated cylinders 
of radii a and b 

Ei2{<j) ^ Ec{(j) - El (<j, a) - El (a, b) . (7) 

The divergences in Ec{c) are cancelled out by those of Ei{(j, a) and Ei{a, b). 
Therefore, to compute Ei2{(t) we can set (f> ^ and a = 0. The contour 



integral reduces to an integral on the imaginary axis. After some steps we 
findHl 

where a = b/a and 
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Eqs. (|HJl and © allow a simple numerical evaluation of E12. The exact energy 
for the two concentric cylinders is the sum of E12 and the Casimir energies 
for single cylinders of radii a and b '^ 

Ec - Ei2 - 0.01356 i\ + ^)Lhc. (10) 

3 Proximity and semiclassical approximations 

The Casimir energy per unit area for parallel plates separated by a distance I 
is given by Eppil) = — fjUF- ^^ ^^^ proximity approximation, the interaction 
energy between two conductors that form a curved gap of variable width z 
can be computed asl^ 

Ei= I Epp{z)d<j , (11) 

where E is one of the two surfaces that define the gap. It is clear that in the 
above approximation the non-parallelism of the area elements is not taken 
into account. Moreover, the result is different if the other surface is chosen to 
perform the calculation. These corrections are expected to be small for low 
curvature and close surfaces. 

For the concentric cylinders, the proximity approximation gives 

As mentioned before, there is an ambiguity in the choice of the effective area 
Aeff. A calculation based on the inner cylinder gives ^e// = ^iraL, while 
for a calculation based on the outer one gives Aeff — 2TrbL. The difference 
is of course harmless in the limit a — > 1. In this limit, the exact result for 
E12 coincides with the proximity approximation. Indeed, using the uniform 
expansion of Bessel functions in Eqs. ^ and ©, one can show that to leading 
order in a — 1, £^12 ~ Ej(b — a). 

For larger values of a one expects the exact result to differ significantly 
from the proximity approximation. We have computed numerically the inter- 
action energy i?i2, and compared it with the proximity approximation, using 



different choices for the effective area. In particular, we parametrized the ef- 
fective area as A^ff = 2TTLaPb^~P, and found that the value p = 1/2 provides 
the best fit for the numerical data. For this value of p, the effective area is 
the geometric mean of the areas of both cylinders. As shown in Fig. 1, the 
proximity approximation with this particular choice for the area reproduces 
the exact results far beyond its expected range of validity: Isl the discrepancy 
for the pressure on the inner cylinder due to the presence of the outer one is 
less than 10% for 1 < a < 4. The choice of the geometric mean is crucial for 
this agreement. 

Let us now discuss briefly the semiclassical approximation for the Casimir 




Figure 1. Dimensionless Casimir interaction energy (upper panel) and pressure (lower 
panel), as a function of a = b/a. In both panels the dashed line corresponds to the exact 
result and the full line to the p = 1/2 proximity result. 



energy, which can be written as 

Ec^ [ l-E posciE) dE , (13) 

Jo ^ 

where Posc{E) is the difference between the spectral density of electromagnetic 
modes in the presence of conductors and the spectral density in vacuum. 
Periodic orbit theory relates posc of a given H arniltonian to the periodic orbits 
in the corresponding classical system. 1 ^ 1 1 I v l To leading order in h, 

P°''{E) = ^ E ^*(^) ^HSt{E)/n + pt) , (14) 

t 

where the sum runs over periodic orbits labeled by t, and St is the classical 
action of the periodic orbit t. For photons of energy _E, St{E) = ELt/c, with 
Lt the length of the periodic orbit. The phase pt is the so-called Maslov 
index. The exponent v and the amplitudes A depend on the type of the 
periodic orbit. 

The periodic orbits in the region between very long concentric cylinders 
are contained in planes perpendicular to the axis of the cylinders. There are 
two different classes of orbits. The type-I orbits do not touch the inner cylin- 
der, and are polygons that may be uniquely labeled by two integers (w,w), 
where v is the number of bounces in the outer cylinder and w the winding 
number around the center. Type-II trajectories do touch both cylinders, and 
are also labeled by {v,w). We have have computed the Casimir interaction 
energy using this semiclassical approximation. 1^ The final result is domi- 
nated by the self-retracing type II periodic orbit (w = 0) and its repetitions. 
Moreover, it coincides with the proximity approximation if the effective area 
is taken as the geometric mean of the areas of both cylinders. 

4 Eccentric cylinders: an experimental proposal 

Up to now we considered concentric cylinders. Obviously, the force between 
them vanishes in this case. Let us now consider two slightly eccentric cylinders 
of radii a and b. We will mainly focus on the particular case a c^ b, since the 
Casimir interaction between cylinders is then stronger. The distance between 
the axis of the cylinders will be denoted by e (see Fig. 2). In order to evaluate 
the Casimir energy for this configuration, we will use the proximity approx- 
imation. This is partially justified by our previous results. The interaction 
energy between cylinders is 

where r{d) is the distance of a point of the external cylinder to the axis of the 
inner one and dA(,ff{9) is the geometric mean of two small adjacent areas on 
both cylinders. From Fig. 2 one can check that r{6) = \/b'^ — e^ cos^ d + esin0 




(a) 




(b) 



Figure 2. Two eccentric cylinders, (a) An inner cylinder of radius a and a hollow cylinder 
of radius b, with the origin of coordinates on the axis of the inner cylinder, and distance e 
between the two axis, (b) The effective area for the application of the proximity approxi- 
mation, as the geometric mean of Si and S2- dA^ff{8) = y/SiS2 



and dAeff — L\/ah + ea sin 9d9. Taking the derivative of the energy with 
respect to e we obtain the force on the outer cyhnder due to the inner one. 
Since we are considering a ~ 6, we will always have e/b <^ 1. Thus to the 
lowest non trivial order we obtain'*'' 
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TT^hcLa 



desine 
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(16) 



24064 J, [i + |sin0_|]4 

where e = e/(6 — a) and Fq = Tr'^hcLa/60{b — a)^. It is worth noting that 
closest surfaces are attracted together, and that the force only vanishes when 
the cylinders are exactly concentric. The equilibrium position is unstable. In 
the particular case in which e ^ 1, the force is linear in the distance between 
the axis of the cylinders Fy « eFo . This corresponds to an inverted harmonic 
oscillator, and explicitly shows the instability . In the opposite case, when 
e -^ 1, the force scales like d~^'^, where d = b — a — e is the minimum distance 
between cylinders. 

We now discuss a possible experimental arrangement that could be used 
to measure the Casimir force between cylinders. '^ We consider the almost 
coaxial configuration e <C 1. The external cylinder could be mounted on a 
resonator of effective mass M and natural frequency ujq. The presence of 
the inner cylinder renormalizes the frequency of the resonator. Assuming 
a small frequency shift we obtain -^ = — 2(b-a)Muj'^ • ^'^^ typical values of 
the parameters it is simple to reach a frequency shift of 0.1%. From an 
experimental point of view, this configuration has some advantages over the 
parallel plates configuration. On the one hand, the alignment procedure is 
easier for the cylinders than for the plates (see below). On the other hand, 
when there is no residual charge in the inner cylinder, the system remains 
neutral and screened by the external one from outer noises, interferences, and 
from residual charges in the outer cylinder. The expected gravitational force is 



obviously null. When the inner cylinder has a residual charge, there will be a 
small potential difference between the cylinders, and the coaxial configuration 
will be electrostatically unstable. To avoid it, one could start the experiment 
by putting both cylinders in contact. Ideally, the charge of the inner cylinder 
would flow to the outer one, minimizing the effect if residual electrostatic 
forces. Alternatively, one could use the electrostatic instability to improve 
the alignment: if a time dependent differential potential is applied between 
cylinders (as in the null experiments to test the exactness of the electrostatic 
inverse square law), parallelism would be optimized when the value of the 
force is minimum. 

To summarize, we have shown that the interaction between cylindrical 
shells is both of theoretical and experimental interest for the study of Casimir 
forces. On the theoretical side, it is a simple but non trivial example for 
testing the validity of different approximations to the Casimir energy. On the 
experimental side, it is a promising geometry for measuring the Casimir force. 
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